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The RENO experiment recently reported the disappearance of reactor electron antineutrinos (i? e ) 
consistent with neutrino oscillations, with a significance of 4.9 standard deviations. The published 
ratio of observed to expected number of antineutrinos in the far detector is R=0.920±0.009(stat.) ± 
0.014(syst.) and corresponds to sin 2 26 13 = 0.113 ± 0.013(stat.) ± 0.019(syst.), using a rate-only 
analysis. In this letter we reanalyze the data and we find a ratio R=0.903 ± 0.01(stat.), leading to 
sin 2 2Si3 = 0.135. Moreover we show that the sin 2 2#i3 measurement still depend of the prompt 
high energy bound beyond 4 MeV, contrarily to the expectation based on neutrino oscillation. 



Introduction 

For ~X-2 km baseline reactor neutrino experiments the 
survival probability of electron antineutrinos (u e ) is 



survival 



1 - sin 2 26»i3 sin 2 (1.27Am 2 L/£0, (1) 



where E is the energy of v e 's in MeV, and L is the 
baseline distance in meters between the reactor and 
detector and Am 2 = (2.32t°; 12 ) x 10" 3 eV 2 0). 

In 2011, first indications of a non-zero 6*13 value have 
been reported by the T2K 0] and MINOS [j] accelerator 
appearance experiments, and by the Double Chooz re- 
actor disappearance experiment [a]. Very recently, the 
Daya Bay experiment reported the most precise mea- 
surement of 0i3 using a rate-only analysis and found 
sin 2 26»i3 = 0.092 ± 0.016(stat) ± 0.005(syst) 0. 
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RENO Oscillation Results 

In a first publication released on April 3 rd 2012 the 
RENO collaboration reported a measurement of the 
neutrino oscillation mixing angle #13, based on the 
observed v e disappearance with a significance of 6.3 
standard deviations [H. They obtained a ratio of 
observed to expected number of v e in the far detector 
of R = 0.922 ± 0.010(stat.) ± 0.008(syst.). Using 
a rate-only analysis, they derived a very precise mea- 
surement of sin 2 26»i3 = 0.103±0.013(stat.)±0.011(syst.). 

On April 8 th 2012, in a second version of their publica- 
tion, the RENO collaboration revised their results [2j. 
The ratio of observed to expected number of £* e 's in the 
far detector has been updated to R=0.920±0.009(stat.)± 
0.014(syst.), though the E prompt energy range being used 
for the computation is not explicitely mentioned. This re- 
vision implies a modification of the measurement of #13 



FIG. 1. Original figure taken from Q|. Upper panel: spectrum 
of the prompt signal observed in the far detector compared 
with the non-oscillation prediction from the £> e spectrum ob- 
served in the near detector. The backgrounds shown in the 
inset are those subtracted to the far spectrum. Errors are 
statistical uncertainties only. Lower panel: ratio of the two 
spectra displayed above. 



leading to sin 2 26*13 = 0.113 ± 0.013(stat.) ± 0.019(syst.). 
It is worth noting that estimate of the total systematic 
uncertainty increased by a factor of 1.7 between the two 
release versions. The collaboration updated the results 
by explicitly including the global normalization in the 
oscillation fit (constrained by the absolute reactor neu- 
trino flux uncertainty of 2.5%) and modifying the lithium 
background determination. The detail characteristics of 
the RENO experiment can be found in [H, Hj] . 
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Expected disappearance signal 

In order to assess the expected disappearance signal at 
the far detector we perform a simulation of the RENO 
experiment. We use the fluxes published in [l| to qualify 
our simulation setup with six nuclear cores and two detec- 
tors, within the percent level. Table Q] shows the RENO 
expected cumulative v e deficit, (5(£ , hi g h), as a function of 
the high prompt energy bound, Shigh- It is given by the 
formula 



lMcV 



P. 



survival 



{E v {E)))dE 



&yS(E)dE 



(2) 



where S(E) is the expected prompt energy spectrum at 
the far detector, taking into account the the six reactor 
distances, powers and fuel compositions. It shows that 
roughly 98% of the v e deficit ('the oscillation signal') is 
expected below 6 MeV (prompt energy). The oscillation 
analysis results must be thus stable by changing the high 
energy bound to values between 6 and 12 MeV. 
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FIG. 2. Weighted averaged of the RENO Far/Near data (bot- 
tom frame of Figure[TJ as a function of the high energy bound 
(statistical errors only). The red points shows the expected 
behavior based on our RENO simulation (sin 2 (26>i 3 )=0.113), 
whereas the black points represent our RENO reanalysis. Our 
measured Far/Near ratio is 20% lower than the RENO pub- 
lished value, R=0.920. Moreover, the data do not match well 
expected evolution for neutrino oscillation. 



TABLE I. Simulation of the RENO cumulative v R deficit in 
the prompt energy range [1;-Ehigh] MeV, assuming Am 2 = 
2.32 x 10~ 3 eV 2 . Most of the oscillation signal is expected 
below 5 MeV. 



Re-analysis of RENO data 

In this section, we present a reanalysis of the RENO 
published data. We use the Far/Near ratio spectrum 
(bottom frame of Figure [T]) . We recompute the mean 
Far/Near ratio over the E prompt = [1-6] MeV energy 
range, and unexpectedly obtain a lower Far/Near ratio 
R=0.903 ± 0.010(stat.), 2% downward shifted with re- 
spect to the RENO published version. We then evaluate 
the impact of this bias on the oscillation result. We simu- 
late the relation between the Far/Near ratio and sin 2 26>i3 
for a rate-only analysis, leading to the expression: 



sin 2 26»i3 



1- < P. 



surviva 



i > 



l-R 



< sin 2 (Am 2 L/4E) > 0.715 



(3) 



valid for -Ehigh >5 MeV. Our reanalysis of the RENO 
data in the E prompt = [1-6] MeV energy range lead to 
sin 2 2^i3 = 0.135. Consequently, the data indicates a 
possible underestimation of sin 2 2#i3 by about 20% (la). 

We now test the consistency of the data points and their 
error bars (bottom frame of Figure [lj . The evolution 
of the Far/Near ratio as a function of the high energy 
bound of the prompt energy range, Ehigh is displayed in 



Figure [U Consequently we study the evolution of the 
sin 2 (26*13) measurements as a function of E] lig h. We re- 
trieve the data from the upper panel of Figure [TJ present- 
ing the far detector spectrum and the expected near de- 
tector non-oscillating spectrum. The upper panel shows 
background subtracted spectra up to 12 MeV and then 
allows us building the Far /Near ratio as a function of the 
prompt energy up to 12 MeV. In the [1-8] MeV energy 
range we cross check our results using the top and bot- 
tom frames of Figured] validating the method of properly 
normalizing both near and far spectra based on published 
distances, live times and efficiencies [l], Q. Results arc 
displayed on Figure [3] where the simulation is displayed 
in red and the data reanalysis in black. The evolution 
on sin 2 (26*13), does not match the expectation were 80% 
of the entire v e deficit is expected below 4 MeV. This 
implies that the measurement of sin 2 (26*13) should not 
significantly change with an higher energy bound con- 
trary to the pathologic behavior of the data leading for 
instance to a modification of sin 2 (26 i i3) value from 0.11 
to 0.135 between 4 and 6 MeV, whereas the v e expected 
deficit varies by less than 2%. 

Surprisingly the rate-only measurement of sin 2 26*13 
still changes when considering energy bound higher 
than 8 MeV, as opposed to expectations (red curve on 
Figure [3]) . This change is clearly not due to neutrino 
oscillation physics. Using £^^=12 MeV one obtains 
sin 2 26*13 = 0.115, consistent with the value published 
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FIG. 3. Simulation and reanalysis of the sin 2 29n constraint 
as a function of the high energy bound of the prompt en- 
ergy range over which the Far/Near ration is computed (for 
sin 2 (2#i3)=0.113, statistical errors only). We use the data 
published in the upper panel of Figure [1] Accounting 99% of 
the neutrino oscillation signal, at Ehi g h=6 MeV, one obtains 
sin 2 2013 = 0.135. But the measured central value sin 2 2#i3 
still decreases towards 0.115 if we (erroneously) include data 
over the whole prompt energy range, [1-12] MeV, contrarily 
to neutrino oscillation expected behavior. 



the collaboration report a deficit of 8.0% in the far detec- 
tor compared to what is expected in the near detector if 
no-oscillation, our analysis, which considers the expected 
neutrino oscillation energy range (E prompt — [1-8] MeV), 
leads to a deficit of 10.0%. Our reanalysis points out to 
a possible bias of the central value by about +20% (Iff), 
leading to a higher sin 2 26*i3 = 0.135. This new best fit 
result is more than 2er off the Daya Bay central value. 
We then study both the Far/Near ratio and the sin 2 26*i3 
measurements as a function of the high energy bound of 
the prompt energy range, Ehigh- The data shows a patho- 
logic behavior, especially beyond 5 MeV, which could be 
interpreted as an underestimation of the relative system- 
atic uncertainty between the near and far detectors or 
as an inaccurate background subtraction. Finally, we at- 
tempt to understand the origin of this sin 2 26*i3 best fit 
discrepancy and we notice that a fallacious inclusion of 
the data between 8 and 12 MeV in the oscillation anal- 
ysis (though no oscillation should occur in this energy 
range) would lead to sin 2 2f?i 3 = 0.115. We kindly ask 
the collaboration to shed more light on their analysis in 
a detailed publication. 



by the RENO collaboration Q. We conclude that the 
deficit above 6 MeV could be interpreted as Far/Near 
detector response differences or as an inaccurate back- 
ground subtraction. 

As a cross-check of the method developed to reanalyze 
the recent RENO reactor antineutrino results [2J, we re- 
compute the Daya Bay Far/Near ratio based on pub- 
lished information appearing in the Figure 5 (bottom 
frame) of @. The detail characteristics of the experi- 
ment can be found in Q. We find an averaged Daya 
Bay Far/Near ratio of R = 0.937 ± O.Oll(stat) in the 
Eprompt — [1-7-8] MeV energy range, which is in excellent 
agreement with the Daya Bay published result. Further- 
more, we study the Far/Near ratio as a function of the 
high energy bound of the prompt energy range, Ehigh- 
The Far/Near ratio is stable when accounting data with 
Ehigh >5 MeV, consistent with neutrino oscillations. 

Conclusion 

In this letter, we reanalyze the reactor neutrino data re- 
cently published by the RENO collaboration [2|. While 
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